Geometric Diagnostics of Complex Patterns: Spiral Defect Chaos 
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Motivated by the observation of spiral patterns in a wide range of physical, chemical, and biological 
systems we present an approach that aims at characterizing quantitatively spiral-like elements in 
complex stripe-like patterns. The approach provides the location of the spiral tip and the size of 
the spiral arms in terms of their arclength and their winding number. In addition, it yields the 
number of pattern components (Betti number), as well as their size and certain aspects of their 
shape. We apply the method to spiral defect chaos in thermally driven Rayleigh-Benard convection 
and find that the arclength of spirals decreases monotonically with decreasing Prandtl number of 
the fluid and increasing heating. By contrast, the winding number of the spirals is non-monotonic 
in the heating. The distribution function for the number of spirals is significantly narrower than 
a Poisson distribution. The distribution function for the winding number shows approximately an 
exponential decay. It depends only weakly on the heating, but strongly on the Prandtl number. 
Large spirals arise only for larger Prandtl numbers (Pr > 1). In this regime the joint distribution for 
the spiral length and the winding number exhibits a three-peak structure indicating the dominance 
of Archimedean spirals of opposite sign and relatively straight sections. For small Prandtl numbers 
the distribution function reveals a large number of small compact pattern components. 



Many systems in nature exhibit complex pat- 
tern that may be stationary or time-dependent, 
possibly in a chaotic fashion. To understand these 
patterns and any transition they might undergo 
it is important to have quantitative measures that 
characterize the relevant properties of the pat- 
terns and their time dependence. Due to the 
multitude of different types of patterns it is not 
to be expected that a single measure would be 
sufficient to capture the qualitatively different as- 
pects of the various patterns. Thus, while quite a 
few different measures have been developed over 
the years, so far no convincing approach is avail- 
able that extracts the characteristic features of 
patterns dominated by spiral-like pattern com- 
ponents. This paper presents a method that al- 
lows to determine quantitatively features like the 
length of a spiral and its winding number. The 
approach is, however, not limited to proper spi- 
rals and yields additional insightful measures. 



I. INTRODUCTION 

Complex spatio-temporal patterns abound in nature. 
Many of them contain also chiral elements like spirals 
or spiral segments. Particularly well-known and beauti- 
ful examples are spiral waves in oscillating and excitable 
chemical systems P, |2, 0, 0, § , chiral structures in micro- 
bial growth patterns [auljand spiral defect chaos in ther- 
mally driven convection 0, Q . Other examples are spi- 
rals arising in vibrationally excited waves on the surface 
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of granular materia l llfl , in calcium waves in oocytes [111 ] 
and hippocampus [l2( , and in the agg regation patterns 
of the amoeba Dictyostelium [l3L Il4| . These patterns 
often evolve chaotically in time and arise in transitions 
from simpler structures and may in turn undergo further 
transitions to other, possibly more complex structures. 
To gain detailed understanding of such spatio-temporally 
chaotic states and their transitions it is important to have 
quantitative measures that extract the relevant features 
of the patterns. 

Reflecting the variety of spatio-temporal structures ob- 
served, a host of different measures and approaches have 
been developed over the years to characterize the vari- 
ous types of spatio-temporal chaos. Global measures for 
the disorder in the pattern can be obtained from the cor- 
relation function and the spectral entropy ^5|. These 
measures are quite general; they do not give any local 
information and do not take any specific features of the 
patterns into account. For patterns that can be con- 
sidered as deformed periodic patterns a set of disorder 
measures has been introduced that is sensitive to the 
curvature and the deformation of the pattern [l6L Il7| . 
In many cases the patterns are locally stripe-like. For 
such patterns approaches have been developed that ex- 
tract the local wavevector Hll. the orientation of the 
stripes [2(|, and their curvature |20| . 

Many complex patterns exhibit striking, isolated ob- 
jects like point defects, e.g. dislocatio ns pH l22 | . discli- 
nations [23 or penta-hepta defects pil |25|. Dislocations 
often take on the form of spirals, as is the case for points 
of vanishin g os cillation amplitude in spatially extended 
oscillations 129, |22| . Similar spirals also arise in excitable 
systems 0, l28j . Spiral-defect chaos in convection simul- 
taneously exhibits dislocations, disclinations, and spirals 
. In systems displaying such point defects it has proven 
useful to investigate the statistics of the defects or spirals 
[Hill IHEl Hi and also of their trajectories pa.l30j. 



2 



Recently, a different type of characterization has been 
developed that focuses on the topology of the patterns 
and extracts the number of contiguous components and 
of their holes. It has been applied to numerical sim- 
ulations of spatio-temporal chaos in excitable reaction- 
diffusion models |3lJ and to extensive experiments on 
spiral defect chaos |3^ |. 

To identify specifically spirals and to measure their 
characteristics only few methods have been used so far. 
In oscillatory systems in which the spirals emanate from 
a defect with vanishing oscillation amplitude the defect 
is relatively easily found and the number of defects can 
be counted [2(|. In [3i| a fast method is discussed to 
find the focus of spirals, independent of their dynamical 
origin, which is based on the evolute of the emitted wave 
pattern. In principle, spirals can also be detected from 
singularities in the local orientation of the pattern 20] or 
the local wavevector [l8j. 

In this paper we discuss an approach that allows to ex- 
tract typical spiral features like the position of the spiral 
tip, the size of the spiral measured in terms of the length 
of the spiral arm as well as its winding number. The 
method is, however, not limited to proper spirals, but 
characterizes arbitrarily shaped segments of stripe-like 
patterns that have a typical wavelength. In the process 
other measures like the number, size, and compactness 
of the components of the patterns are obtained, as well. 

We discuss our diagnostic method at length in SecllTl 
In Sec lIIII we apply it to spiral defect chaos in thermal 
convection and present the dependence of various mea- 
sures on the heating and on the Prandtl number of the 
fluid. We show that the distribution function for the 
number of spirals in the pattern is narrower than a Pois- 
son distribution. The distribution of the winding number 
of the spirals is reasonably approximated by an expo- 
nential as had been suggested by experiments j34[- We 
find that while the decay rate depends only weakly on 
the heating it does depend strongly on the Prandtl num- 
ber. Correlating the contour length and the area of the 
components we find the signature of target patterns for 
Pr = 1.5, but not for Pr = 0.3. For Pr = 1.5 the joint 
distribution for the spiral length and the winding number 
reveals that Archimedean spirals and relatively straight 
components dominate the patterns, while for Pr = 0.3 
no such distinct groups are found. Instead, in this regime 
the patterns are characterized by a large number of small 
compact components. In Sec llVI we discuss our conclu- 
sions. 



II. METHOD 

Our main goal is to identify spiral-like structures in a 
sequence of patterns and to characterize them quantita- 
tively. In the process we also obtain some other measures 
of the patterns. Here we describe the procedure in detail. 
It is based on a set of contour lines giving the pattern. 
Spiral arms are defined via the tip at the point of highest 



curvature and the merging of the arm with other parts 
of the pattern at a vertex (cf. Fig^ below). 

We start from a sequence of snapshots that character- 
ize some continuous quantity of the system. In our exam- 
ple of convection it is the temperature field in the mid- 
plane of the convection cell. While in many pattern anal- 
yses such a continuous representation is first reduced to 
a binary representation, we use the continuous represen- 
tation to obtain the set of contour lines corresponding to 
the level midway between the minimal and the maximal 
value of the field taken across all frames of the sequence. 
The resulting contour lines retain more spatial resolu- 
tion than the interface of a binary pattern. Our simula- 
tions of the convection system employ periodic boundary 
conditions. To reflect this periodicity, we appropriately 
connect contours that reach the boundaries. Most con- 
tours become closed in this process. In many patterns 
there are, however, also a few open contours; they wrap 
around the whole system and may, for instance, corre- 
spond to a straight convection roll that spans the system 
and connects with itself across the periodic boundaries. 

Having obtained all closed contours we already are in 
the position to obtain a few simple quantitative charac- 
terizations of the pattern. The number of closed con- 
tours provides a topological measure of the pattern by 
giving its number of components. In many patterns it is 
useful to distinguish between 'white' and 'black' compo- 
nents, e.g. contours inside of which the temperature is 
higher (lower) than the level of the contour. Focussing, 
for instance, on the pattern made up of the white com- 
ponents, the number of white components gives then the 
Betti number of order 1, while in our case of periodic 
boundary conditions the number of black components is 
related to the Betti number of order 2, which counts the 
number of holes in a pattern j^, [3^. It is often use- 
ful to go beyond the topological characterization of the 
patterns and also measure metric aspects like the size of 
components. We therefore keep track also of the length 
V of the perimeter of each component, i.e. the length of 
the enclosing contour, and of its area A. Comparing the 
length and area of a given closed contour gives a measure 
for the compactness of the component, which we define 
as 



With this normalization a circle has C — 1, while small 
values of C correspond to filamentary structures. 

To extract spirals from the contour lines we start by 
finding the tips of the spirals, which we identify with 
the points of highest positive (convex) curvature of the 
contour. In disordered patterns and, in particular, in 
noisy patterns points with high curvature may also ap- 
pear away from the tips. To reduce the number of these 
spurious tips we first smooth the contour. To do so we 
'diffusively' move the polygon points (xi 7 yi), i = 1 . . . N, 
which define the contour, using the mapping 

Xi — > x l + 5(x i+ i - 2x l + Xi-i) (2) 
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FIG. 1: a) Example of a convection pattern for Pr = 1.0 
and e — 1.0 exhibiting spiral defect chaos, b) Example of a 
component in the left-bottom quadrant of the pattern in a). 
Original contour (green open circles), smoothed contour (red 
open circles), local maxima and minima of the curvature (red 
and blue solid squares, respectively), inflection points (green 
solid squares), shortest distance across the components (blue 
line). Here the axes are in terms of the grid coordinates rather 
than physical coordinates. 

and analogously for yi. To avoid overshoots we use a 
small value of 5, typically 8 = 0.2, and perform the 
smoothing iteratively (typically in 50 steps). While the 
smoothing strongly reduces localized bulges along overall 
straight sections of the contour and therefore reduces the 
number of spurious tips, the curvature at a spiral tip is 
not affected much. Note that small components shrink 
substantially by the smoothing. In this procedure we 
are, however, mostly interested in large components and 
their reduction in size is quite small. This is illustrated 
in Fig^D which depicts a white component in the bottom 
left quadrant of the pattern shown of Fig^Ji before (green 
open points) and after (red open circles) the smoothing. 
A further benefit of the smoothing is that it distributes 
the grid points {xi,yi) more equally along the contours. 

To obtain the local curvature of the contour we 
parametrize for each point (xi, yi) of the contour a small 
section of the contour consisting of 5 to 11 points around 
(xi, yi) by the arclength s and perform a best fit with two 
quadratic polynomials (x(s),y(s)). The local curvature 
is then given by 

_ x'x"-y'y» 
y/x' 2 + y' 2 

where the primes denote first and second deriva- 
tives of x(s) and y(s), respectively, at the midpoint 
(x(s m ), y(s m )) of the arc approximating the contour over 
this local range. 

We then identify the spiral tips contained in this con- 
tour iteratively by locating the point with maximal pos- 
itive curvature. To measure the length of this spiral arm 
we follow both branches of the contour originating from 
the tip until the first vertex is reached (cf . Fig^) . Ver- 
tices are identified by the width of the arm exceeding a 
threshold w max . Typically we take a value of w max = 10, 
which corresponds to about one wavelength of the pat- 
tern. To measure the width we find for each point on one 



branch emanating from the spiral tip the closest point on 
the other. In general this is a non-local problem. It be- 
comes much simpler by starting from the tip (xi tiv , y% u ) 

and finding for each consecutive point P- = {pyp ,y. ), 

j < Hip, on the 'lower' branch the point Pj. u \P^) on the 

'upper' branch that is closest to P- and that is within 
a certain range i t i P < k' < k < k' + Ak of the previously 

identified point P fe ( , u) (Pf\) on the upper branch. We typ- 
ically take Ak — 5. In Fig^> the closest points identified 
in this manner are connected by blue lines. Near the tip 
this procedure sometimes leads to connections that lie 
outside rather than inside the contour. To aid avoid- 
ing this problem, we require that within a core region 
near the tip the point on the upper branch strictly ad- 
vance, k > k' . Typically we take the core region to be 
10 points. When reaching a vertex, typically the points 
on the lower branch still advance while those identified 
as closest on the upper branch do not (e.g. near x = 115 
y = 197 in Fig|TJ)). We then identify the end of the 

spiral arm with the last pair of points (P- ,Pj. u (Pj )) 
for which the point on the upper branch still advanced, 

P^ u \Pj [S> ) 7^ Pfr (Pj l 2i)- For a segment or spiral arm 
thus identified we take its length S to be the mean of the 
arclength along its two branches and define its winding 
number via W = + 9^)/2n, where <?« and 

are the angles by which the two branches of the contour 
are rotated from the tip to the vertex, respectively. 

Having found the first tip and the associated spiral arm 
we then identify the second spiral tip with the point of 
maximal positive curvature among the points on the con- 
tour not included in this first arm and correspondingly 
for further arms. When identifying the subsequent spiral 
arms we ensure that they do not overlap with previously 
identified arms. In many contours there are portions that 
are not associated with a spiral tip. Nevertheless, even 
in these parts of the contour the algorithm may identify 
segments that are connected with local positive maxima 
of the curvature. Fig^D shows such a segment in the 
lower right. It is associated with a black inclusion. The 
size of the segments generated in this way is of the order 
of w max - To exclude such segments, we omit all seg- 
ments with an arclength smaller than 2w max in the final 
analysis of the results. 

Many contours do not include a vertex. In those cases 
we associate the contour with two segments (spiral arms) 
by cutting the contour at the point P^ (and P^ U \Pj^)) 
at which the absolute value of the curvature outside the 
core regions of the two tips is minimal. 

The spirals of spiral defect chaos in convection are not 
topologically stable objects; their winding number varies 
continuously and is not conserved. In other systems spi- 
rals are topologically stable and can only be created in 
pairs with opposite topological charge. Well-studied sys- 
tems of this type are oscillations in extended media in 
which the spiral cores represent singularities of the phase 
of oscillation at which the oscillation magnitude vanishes 
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[26(. In particular in such cases it is of interest to count 
the number of spirals with positive and negative charge 
independent of their size. Employing the algorithm dis- 
cussed here, the topological charge of the spirals can be 
related to the location of inflection points of the con- 
tour in the vicinity of the spiral tip. Even in quite small 
spirals the curvature changes sign near the tip and the 
side on which this inflection point occurs depends on the 
topological charge of the spiral. Note, however, that even 
segments that have a noticeable winding number may not 
be classified as spirals according to this criterion if their 
tip is part of a relatively straight section. An example of 
such a segment is seen in Fig^i where the segment with 
winding number W ~ 0.25 is quite straight near its tip 
and consequently the first inflection point (marked with 
a green solid circle) is relatively far from the tip. 



III. GEOMETRIC ANALYSIS OF SPIRAL 
DEFECT CHAOS 

We use the approach discussed in Sec HTIt o analyze 
spiral defect chaos in Boussinesq convection 0, || , which 
is characterized by two system parameters: the Prandtl 
number Pr of the fluid, which measures the ratio between 
its viscosity and its heat conductivity, and the heating, 
which is measured by the dimensionless Rayleigh number 
R. The latter is often expressed in terms of the scaled 
distance e from the onset of convection at R c = 1708, 
€ = {R-R C )/R C . 

Our results are based on direct numerical simulations 
of the Navier-Stokes equations employing a pseudospec- 
tral code developed by W. Pesch and co-workers |35| . 
Details of the equations and the code can be found in 
[3ri |3tJ ■ The code is based on Fourier modes in the hori- 
zontal direction and appropriate combinations of trigono- 
metric and Chandrasekhar functions that satisfy the top 
and bottom boundary conditions in the vertical direc- 
tions |33,|39|. All runs are performed with 256 x 256 hor- 
izontal Fourier modes for a system size of L — 32 • 2i:/q c , 
where q c = 3.098 is the critical wavenumber, and 2 ver- 
tical modes. Previous simulations have shown that this 
spatial resolution is sufficient to capture the spiral defect 
chaos at the very least semi-quantitatively [3Jf, 36] . Since 
the data analysis requires relatively long runs, going to a 
significantly higher spatial resolution is beyond our cur- 
rent computational means. To solve for the time depen- 
dence the code uses a fully implicit scheme for the linear 
terms, whereas the nonlinear parts are treated explic- 
itly using a second-order Adams-Bashforth method. The 
time step is typically taken to be i^/500, where t v is the 
vertical diffusion time. In our analysis we use snapshots 
of the temperature field in the mid-plane of the convec- 
tion cell that are typically taken every 2t v . We focus on 
three different values of e and Pr each. The number of 
snapshots evaluated range from 500 for Pr = 0.3 to 2000 
for Pr = 1.5. Examples of patterns at e = 1.0 are shown 
in Figflji for Pr = 1.0 and in Fig|2K,b for Pr = 0.3 and 




FIG. 

Pr = 



2: Snapshots of convection patterns at e = 1 for a) 
0.3 and b) Pr = 1.5. System size L = 64.9. 
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0.5 1 1.5 

Red. Rayleigh Number e 

FIG. 3: Mean number of closed contours as a function of e for 
different values of the Prandtl number Pr. 



Pr = 1.5. 

We first discuss the results for various mean quanti- 
ties. FigOshows the dependence of the mean number Af 
of closed contours on e for three different Prandtl num- 
bers. More precisely, Af gives the mean (Afb + Af w )/2 of 
'black' and 'white' components. To assess the variations 
in Af over the course of the simulation Fig|3| shows in 
small symbols the mean values obtained from averaging 
only over the initial, middle, and final third of the run. 
In most cases the variations are so small that the addi- 
tional symbols are barely visible. The number of closed 
contours increases with increasing e and decreasing Pr. 
This is consistent with the fact that the driving force of 
spiral defect chaos is a large-scale flow, which is driven 
by the curvature of the convection rolls and which in- 
creases with increasing e and decreasing Pr. Such an 
e-dependence has also been found in recent experiments 
for Pr ~ 1 |32| . There the Betti numbers of order 1 and 2 
have been measured, which give the number of contiguous 
components and of holes in the pattern, respectively |3lj |. 
The mean length of the contours exhibits the opposite be- 
havior (data not shown) . In Fig0] the analogous trend 
of increasingly finer structures with increasing strength 
of the large-scale flow is depicted in terms of the mean 
arclength of the segments (spirals) of the components. 

The arclength results for Pr = 1.5 show a large scat- 
ter at e = 0.7. This reflects the fact that for this system 
size e = 0.7 is close to the onset of spiral defect chaos 
[icj . FigO shows that in this regime the pattern does 
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0.5 1 1.5 

Red. Rayleigh Number e 

FIG. 4: Mean arclength 5 of the segments of the components 
as a function of e for different values of the Prandtl number 
Pr. 



not exhibit much of a spiral character. Our computa- 
tional resources do not allow to establish whether this 
low chaotic activity persists indefinitely. Alternatively, 
the pattern could eventually freeze in a disordered pat- 
tern. Experimentally, it has been found that near the 
transition to spiral defect chaos the system can exhibit 
very intermittent behavior in which long periods of al- 
most ordered patterns alternate with states exhibiting a 
few spirals |4f| . 

To extract the spiral character of the chaotic patterns 
we measure their winding number. Since the system 
is chirally symmetric the mean winding number is very 
close to in all cases. The spiral character is therefore 
better captured by the standard deviation of the wind- 
ing number. Its dependence on e and Pr is given in 
FigEl Two trends can be discerned. For the strongly 
chaotic states, i.e. for low Prandtl number and large e, 
the standard deviation of the winding number decreases 
with increasing e and decreasing Pr, similar to the de- 
pendence of the arclength. This reflects the fact that 
shorter spiral segments have a smaller winding number. 
For low chaotic activity, however, the winding number 
decreases with decreasing e, even though the arclength 
continues to increase with decreasing e. This is very 
striking for Pr — 1.5 and e = 0.7 and consistent with 
the visual appearance of the pattern shown in Fig [5] But 
even for Pr — 1.0 the spiral character of the state is non- 
monotonic in e. For these parameter values no tendency 
for the pattern to become ordered was apparent visu- 
ally. Thus, the standard deviation of the winding num- 
ber may be suitable as a sensitive quantitative measure 
for the transition to spiral defect chaos, complementing 
the previously used mean number of spirals [4(j ■ 

A natural conjecture is that the dependence of the 
various measures on e and Pr would be simpler if they 
were considered as functions of the distance from the on- 
set of spiral defect chaos, e — £sDc(P r )- To test this 
in detail would require a quantitative measurement of 
£SDc{P r )- This is beyond our computational capabil- 
ities. However, while both the mean number of com- 
ponents and the arclength take on comparable values 
for {Pr = 0.3, e = 0.4), (Pr = 1.0, e = 1.0), and 
(Pr = 1.5,e = 1.4), this is not the case for the stan- 




FIG. 5: Snapshot of a typical convection pattern for Pr = 1.5 
and e = 0.7. 
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FIG. 6: Standard deviation of the winding number as a func- 
tion of e for different values of the Prandtl number Pr. 



dard deviation of the winding number, which reaches 
significantly higher overall values for Pr = 1.5 than for 
Pr = 0.3. Thus, the values at different Prandtl numbers 
cannot be related to each other by a simple Pr-dependent 
shift and rescaling of e, indicating a somewhat more com- 
plex scenario (see also Figs ll2ll31 below 1 ). 

In systems in which spirals constitute topologically 
conserved defects it has proven quite instructive to mea- 
sure the distribution function for the number of spirals in 
a given snapshot 0, |2^, [2|| [2^, [2j| . In various, but not 
all [2gl |29| systems the distribution function was found 
to be well fit by a squared Poisson distribution, which is 
consistent with a simple over-all behavior of the spirals: 
they are created pairwise with a fixed probability, diffuse 
randomly, and annihilate each other pairwise upon colli- 
sions. The spirals in spiral defect chaos are not topolog- 
ically conserved and may exhibit different statistics. In- 
deed, in a substantial effort the distribution function was 
measured previously by visual inspection of many experi- 
mentally obtained snapshots filj . There it was found for 
Pr = 0.96 in a system with aspect ratio T = 52 that the 
probability distribution function was consistent with a 
Poisson distribution for smaller values of e (e = 0.72 and 
e = 0.80), while for e = 0.96 the distribution was more 
sharply peaked than a Poisson distribution. Of necessity, 
the manual counting limited the number of snapshots 
that could be analyzed and the results have a somewhat 
larger margin of uncertainty. In addition, the distribu- 
tion was determined only for spiral arms with a minimal 



6 




FIG. 8: Rescaled distribution function for the number of spi- 
rals for the three different values for the threshold W mm • The 
curves are obtained from those of Fig[7|by shifting and rescal- 
ing by the mean of the distribution. 



winding number of W m in = 1- To address the question 
of the distribution function more precisely and its depen- 
dence on the selection criterion for the spirals we show in 
Fig0the distribution function for three different values 
of the threshold, W min = \, W min = |, and W m , n = 1 
(Pr = 1.0 and e = 1.0). Our results confirm with higher 
accuracy the experimental finding that the distribution 
function is narrower than a Poisson distribution. More- 
over, shifting the distribution functions by their means 
and rescaling them by the square root of their means, 
which would correspond to the standard deviation for a 
Poisson distribution, results in the distributions shown 
in Fig [H] (note the logarithmic scaling). Within the sta- 
tistical error the distribution functions collapse indepen- 
dent of the threshold W m in and are significantly more 
sharply peaked than the Poisson distribution. So far no 
convincing model has been put forward that provides an 
understanding of the form of these distribution functions. 

Further insight into the patterns and their dependence 
on the system parameters can be gained by considering 
correlations between various properties. An instructive 
example is the joint distribution function for the log- 
arithm of the compactness of the closed contours and 
the logarithm of their contour length, which is shown in 
Fig|5] We find that most components have a compact- 
ness C ~ 1 and small contour length. For larger contour 
lengths a relatively straight ridge is found in the distri- 
bution function. It reflects the fact that the stripes have 
a typical width A/ 2. Thus, for large values of V the com- 




FIG. 9: Joint distribution function for the logarithm of the 
compactness, InC, and the logarithm of the contour length, 
InP, for e = 1.0. a) Pr = 1.5, b) Pr = 1.0, c) Pr = 0.3. 

pactness C can be approximated by 



A few large components deviate from this relationship 
and have a larger compactness; no component was found 
that has a smaller compactness. This asymmetry arises 
because the wavelength of the pattern cannot become 
smaller than a certain minimal value. For large Prandtl 
number an additional structure is apparent: for values of 
the compactness C < 1 the distribution function shows 
three quite distinct peaks in its dependence on the con- 
tour length. We associate these peaks with the appear- 
ance of target patterns, which are comprised of a discrete 
set of quite compact contour lines. When the Prandtl 
number is decreased these peaks become less pronounced 
(cf. Pr — 1.0 in Fig^) and disappear for Pr = 0.3 
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FIG. 10: Correlation between arclength and winding number 
for e = 1.0. Each dot represents one segment of a component. 
The line is a fit to an Archimedean spiral (cf. 10). 



(FigEt). Instead the main peak becomes quite broad and 
extends towards small contour lengths indicating the ap- 
pearance of relatively small compact components. Com- 
paring the distribution functions for different Pr it is also 
apparent that the relative frequency of large components 
decreases noticeably with decreasing Prandtl number. 

Direct information about the spiral character of the 
patterns is obtained by considering the correlation be- 
tween the winding number and the arclength of the seg- 
ments. This is shown in Fig llOl for three Prandtl num- 
bers at e = 1.0. Since the wavelength of the stripes is 
essentially fixed, the winding number W of any spiral is 
expected to be limited by that of an Archimedean spiral 
with the same arclength S. A simple computation gives 

S = irWay* 'Att 2 W 2 + 1- -ahxa 



+ ialn (2naW + a\/ ' An 2 W 2 + lj 



(5) 



As FigfTUI shows. eq.(|5j) fits the extremal values of W 
quite well for all three Prandtl numbers. The fitted value 
of a — 1.8 corresponds to a wavenumber q — 2.2 for the 
rolls making up the spiral, which is consistent with the 
peak of the azimuthally averaged wavevector spectrum 
measured experimentally Q . A cut through the distribu- 
tion function covering the range 10 < 5 < 12.5 is shown 
in Fig^2f° r three values of Pr. For Pr = 1.5 it has three 
peaks. The outer two correspond to near-Archimedean 
spirals while the center one to quite straight structures. 
As the Prandtl number is decreased the triple-peak struc- 
ture disappears and no signature of the Archimedean spi- 
ral is left. As the arclength S is increased the distribution 
function drops off very rapidly for Pr — 0.3 and, as men- 
tioned before, spirals with winding numbers well above 1 
are found only for large Prandtl numbers. 

In experiments using CO2 with Pr = 0.98 it has been 
found that for e = 0.72 and e = 0.84 the distribution 
function for the winding number decays approximately 
exponentially |4l]][5(|. Due to the manual counting of 
the spirals only a limited amount of data was available. 
We have extracted this distribution function also from 
our data. In Fig^| we show the distribution for the 
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FIG. 11: Distribution function of the winding number for 
arclengths S in the range 10 < 5 < 12.5 (cf. FigEUJ. 
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FIG. 12: Distribution function of the winding number for 
Pr = 1.0 and different values of e. 



magnitude of the winding number for three values of e 
for Pr — 1.0. In particular for e = 0.7 quite convinc- 
ing exponential behavior is found. With increasing e the 
function decays, however, somewhat faster at large wind- 
ing numbers, casting some doubt on the exponential be- 
havior. Overall, there is, however, little change with e 
and fits to exponentials would lead to similar decay con- 
stants. By contrast, when considering different values of 
the Prandtl numbers quite different results are obtained 
as shown in Fig ll3l where we keep e fixed at e = 1.0. 
With increasing Prandtl number the decay of the distri- 
bution function becomes significantly slower and spirals 
with larger winding numbers become much more likely, as 
might have been surmised from Fig llOl Thus, as became 
apparent already in the discussion of Figs 1^14161 the state 
of spiral defect chaos depends on the Prandtl number and 
e separately and not only on the distance e — esDc(Pf) 
from its onset. 

Since the spirals break the chiral symmetry a particu- 
larly interesting question is how they respond to an exter- 
nal breaking of the chiral symmetry as it arises in a rotat- 
ing convection system. This question has been addressed 
previously in experiments using CO2 at Pr = 0.98. 
There the number of clock-wise and counter-clockwise 
spirals was counted by hand and it was found that the 
rotation introduces an imbalance between the two popu- 
lations 01 ■ I n terms of the winding number we expect a 
smooth shift of the distribution function with increasing 
rotation rate f2, as is indeed borne out in Fig^] The re- 
sulting dependence of the mean winding number on the 
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FIG. 13: Distribution function of the winding number for 
e — 1.0 and different values of Pr. 
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FIG. 16: Distribution function of the arclength of the spiral 
segments for different rotation rates (Pr = 1.0 and e — 1.0). 
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FIG. 14: Breaking the chiral symmetry. Distribution func- 
tion of the winding number for different rotation rates of the 
system (Pr — 1.0 and e = 1.0). 



rotation rate is shown in Fie 1151 

As shown in Fig ll5l with increasing |0| the mean wind- 
ing number decreases sublinearly. This is connected 
with the fact that the segments of the pattern become 
shorter with increasing rotation rate in anticipation of the 
Kuppers-Lortz instability This decrease is shown 

in Fig ll6l which presents the distribution function for 
the arclength for different rotation rates. Note, that for 
Pr ~ 1 the onset of the Kiippers-Lortz instability is not 
reached until O ^ 14 9J. If instead of the mean winding 
number the mean topological charge of the patterns is 
plotted an essentially linear behavior is found (data not 
shown). 
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FIG. 15: Dependence of the mean winding number on the 
rotation rate for Pr = 1.0 and e = 1.0. Small symbols indicate 
the means over the initial, middle, and final part of the run, 
respectively. 





z=-0.25 


z=0.25 


z=-0.125 


z=0.125 


z=0 


# Contours ('white') 


1.40 


0.733 


1.14 


0.876 


7.15 


# Contours ('black') 


0.738 


1.32 


0.896 


1.13 


7.28 


# Contours 


1.065 


1.031 


1.017 


1.002 


7.21 


# Spirals 


1.09 


1.11 


1.02 


1.06 


1.22 


Arclength 


0.996 


1.03 


0.993 


1.02 


48.07 


SD Winding Number 


1.01 


1.05 


1.001 


1.02 


0.512 



TABLE I: Dependence of the mean number of closed contours 
(positive, negative, and average of the two (cf. FigEJ), mean 
number of spirals, the arclength, and the standard deviation 
of the winding number on the vertical position of the tem- 
perature field. The first four columns of the table give the 
ratio between the quantity at the indicated z-level and at the 
midplane (z — 0). The last column gives the value at the 
midplane (Pr = 0.80, e = 1.4). 



So far we have focussed in our analysis on the tem- 
perature field in the midplane z = of the convective 
layer. While it seems reasonable to assume that the pat- 
tern in the midplane is representative for the convective 
patterns, we test this assumption explicitly by measur- 
ing the mean number of closed contours ('white' compo- 
nents, 'black' components, and the mean of the two), the 
mean number of spirals with W > 1, the mean value of 
the arclength, and the standard deviation of the winding 
number also at different heights z. For these quantities 
Table[|]gives the ratio between the value at the indicated 
z-value and the value at the midplane z = 0, which is 
given in the last column. In these runs we use a reduced 
system size of L — 16 • 2ir/q c retaining 128 x 128 Fourier 
modes. The Prandtl number is Pr = 0.8. 

Considering the number of contours enclosing 'white' 
components, one notices a strong dependence on the 
height z at which the temperature field is measured. For 
instance, while at z — —0.25 the mean number of white 
components is 40% higher than at z = 0, it is almost 30% 
lower at z — +0.25. However, if one averages over the 
white and the black components this variation is reduced 
to 6.5% or less. The z-dependence of the other quanti- 
ties is similarly low. Since the experimental visualization 
techniques measure some vertical average of the temper- 
ature field one can expect that the difference between 
experimentally determined quantities and the computa- 
tionally obtained ones remains below 10%. In the case 
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FIG. 17: Dependence of the distribution function of the ab- 
solute value of the winding number on the vertical position 
z at which the temperature field is measured (Pr = 0.8 and 
e = 1.4). 



of the arclength and the standard deviation of the wind- 
ing number the difference is expected to be even smaller. 
For some quantities like the average number of closed 
contours the variation with Pr and e over the range con- 
sidered is noticeably larger than these differences. On the 
other hand, some quantities like the standard deviation 
of the winding number vary with e and Pr only by 10% 
or 20%. For these quantities the direct comparison be- 
tween the computation and experimentally determined 
values may be not quite as conclusive. However, even 
if different experimental or computational measurement 
protocols should lead to results that differ quantitatively 
from each other, we expect that even less marked trends 
like the non-monotonicity of the standard deviation of the 
winding number (cf. Fig|5J| should be readily resolvable 
as long as the protocol is kept the same when changing 
e or Pr. 

As a further test of the significance of the z-dependence 
Fig |17l shows the distribution function for the winding 
number at z = 0, z = —0.125, and z = —0.25. Within 
the statistical error these distributions are indistinguish- 
able. We conclude therefore that the measures discussed 
in this paper will be meaningful and will provide valuable 
insight into these complex spiral patterns even though 
the patterns themselves - and with them the measures - 
vary somewhat across the layer. 



IV. CONCLUSION 

Motivated by the fascinating convection patterns found 
in spiral defect chaos we have introduced an approach 
to characterize such patterns based on various features 
of the spirals. In the process we have introduced other 
measures that are not specific to spiral-like structures. 
We have used these various measures to characterize spi- 
ral defect chaos as a function of the heating e and the 
Prandtl number Pr. 

Since the onset esDC of spiral defect chaos depends 
on the Prandtl number one could surmise that the prop- 
erties of the resulting spiral defect chaos state depend 
mostly on the distance e — esDc(P^) from that thresh- 



old. Comparing the e-dependence of the mean number 
J\f of components of the pattern and of the arclength S 
of its segments with the e-dependence of the standard 
deviation of the winding number of the spirals we found 
evidence that this simple relationship is, however, not 
satisfied. If any scaling relationship should hold at all, 
it would have to involve also rescaling of the dependent 
variables. Further support for this conclusion comes from 
our finding that the distribution function of the winding 
number is barely affected by changes in e, but it depends 
significantly on the Prandtl number. 

We find that while the arclength of pattern segments 
increases monotonically as the onset of spiral defect chaos 
is approached this is not the case for the winding number 
of those segments. Thus, further away from threshold the 
winding number increases with the arclength as e is de- 
creased, indicating a lengthening of spirals without much 
change in their structure. However, as the threshold is 
approached the segments still become longer but their 
winding number starts to decrease, i.e. the segments be- 
come straightcr. Possibly, the standard deviation of the 
winding number could therefore serve as an early indica- 
tor for the transition to the stationary state. 

With improved accuracy we confirm that the distri- 
bution function for the number of spirals is narrower 
than a Poisson distribution Moreover, using dif- 

ferent thresholds for the classification of the spirals, we 
find that after rescaling these different distribution func- 
tions collapse into a single curve. So far no theoretical 
explanation for this distribution function is available. 

Our automated analysis also allows higher accuracy for 
the measurement of the distribution function of the wind- 
ing number. In previous experiments a roughly exponen- 
tial behavior had been found |4l|. Overall, we confirm 
this behavior, but at large winding numbers the distri- 
bution functions typically decay somewhat more rapidly. 
Interestingly, as mentioned above, the decay depends 
very strongly on the Prandtl number, but very little on 
e. 

By correlating different measures we presented a more 
detailed characterization of the patterns. For instance, 
correlating the contour length of a component and its 
compactness, which is related to its area, it is quite ev- 
ident that for large Prandtl numbers the patterns have 
many filamentary components and a sizeable number of 
target patterns. For small Prandtl numbers the signature 
of the long filaments and of the target patterns disap- 
pears and instead a much broader spectrum of smaller, 
more compact components arises. In a separate inves- 
tigation |45j we find that the number of small compact 
components strongly increases when the fluid properties 
depend significantly on the temperature, i.e. when non- 
Boussinesq effects become important. They introduce a 
resonant triad interaction between the stripe modes that 
enhances the tendency towards hexagonal (cellular) pat- 
terns. Of course, the breaking of the up-down symmetry 
by the non-Boussinsesq effects leads also to differences 
between closed contours that enclose up- flow rather than 
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down-flow regions ('black' and 'white' components). 

Correlating the arclcngth with the winding number al- 
lows the extraction of the spirals. For large Prandtl num- 
bers not only the very long components but also medium- 
sized components tend to be close to Archimedean spi- 
rals. For small Prandtl numbers, however, even among 
the longer components most have winding numbers 
that are much smaller than the winding numbers of 
Archimedean spirals with the same arclength. 

A strength of the measures employed in this paper is 
that they go beyond purely local properties of the con- 
tour lines and identify larger components making up the 
patterns. It should be noted, however, that local changes 
that reconnect stripes can have large effects on proper- 
ties like the arclength and the winding number of the 
affected pattern components. Moreover, in a given snap- 
shot there are only few large components. As a conse- 
quence of both aspects, these measures fluctuate more 
strongly than purely local measures (e.g. [16|), neces- 
sitating a larger number of snapshots. In particular in 
experiments this should pose no problem. 

Since the approach discussed here identifies individ- 
ual spirals one would expect that it would be suitable 
to follow the evolution of individual spirals dynamically. 
This would allow insight into the mechanisms that set 
the mean size of the spirals and may help understand 
the key elements that maintain the chaotic state. It has 
been discussed that dislocation motion and wavenumbcr 
selection play an important role |46( . Visual inspection of 
the numerical results indicates that spirals shrink largely 
through the compression of their outerparts which then 
undergo a skewed- varicose instability [!j . 

We expect that the measures introduced here can be 
used beneficially to investigate whether there are any 
transitions between different states of spiral defect chaos 



as e (or Pr) is changed. Such an investigation would 
require, however, a much finer scan of the parameter val- 
ues, which is well beyond the reach of a computational 
approach. Recently, such transitions have been reported 
based on experimental investigations using the number 
of contours (Betti numbers) as the sole measure • We 
expect that our measures would provide detailed insight 
into the nature of these transitions. At larger e a transi- 
tion from spiral defect chaos to a chaotic state dominated 
by targets has been observed as the Prandtl number is 
increased beyond Pr = 3.5 ^ would be inter- 

esting to see whether the additional measures introduced 
here shed further light on these aspects of spiral defect 
chaos. 

In conclusion, we have introduced a new set of mea- 
sures for the analysis of complex patterns, focusing on 
patterns that contain spiral-like structures. Our results 
for spiral defect chaos suggest that these measures should 
provide also insight into other complex states that exhibit 
chiral structures. Examples include meande ring spirals 
and defect chaos in chemical systems 0, 0, 0, 128] . spirals 
in verticallyvibrated material |10J. patterns of bacterial 
coloni es Vj Q , calcium waves [ll|, [I3 , spiral waves in the 
heart |49j| and in aggregating amoebae 0, Q . 
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